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Abstract 
When the cosmic microwave photons scatter electrons just prior to the decoupling of 
matter and radiation, magnetic fields do contribute to the Stokes matrix as well as to the 
scalar, vector and tensor components of the transport equations for the brightness pertur- 
bations. The magnetized electron-photon scattering is hereby discussed in general terms by 
^ including, for the first time, the contribution of magnetic fields with arbitrary direction and 

in the presence of the scalar, vector and tensor modes of the geometry. The propagation of 
relic vectors and relic gravitons is discussed for a varying magnetic field orientation and for 
different photon directions. The source terms of the transport equations in the presence of 
the relativistic fluctuations of the geometry are also explicitly averaged over the magnetic 
field orientations and the problem of a consistent account of the small-scale and large-scale 
magnetic field is briefly outlined. 
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1 Formulation of the problem 

The last electron-photon scattering is customarily discussed without the additional complica- 
tion of a magnetic field. In numerical codes as well as in analytical estimates, the collisional 
contributions are evaluated as if electrons and ions were free right before last scattering pQ. 
The relativistic fluctuations of the geometry are included in the classic transport problem 
PQ either by using specific gauges [2] or with fully gauge-invariant methods. The resulting 
equations including both the source terms coming from electron-photon scattering and the 
relativistic fluctuations of the geometry form the set of transport equations which can be 
solved within various approaches either by truncating the system at a specific (maximal) 
multipole j2] or by using the integration along the line of sight [3]. One of the consequences 
of the consistent solution of the system of transport equations are estimates of the tem- 
perature and polarization inhomogeneities of the Cosmic Microwave Background (CMB in 
what follows). The recent WMAP 7 data [51 EH El El OH] are able to constrain the vanilla 
ACDM scenario (where A stands for the dark-energy component and CDM for the cold dark 
matter component). In the near future the ACDM scenaridj will be tested not only in its 
minimal version but in its non-minimal extensions ranging from the addition of a stochastic 
background of relic tensor modes of the geometry to large-scale magnetic fields [HI [12]. 

In recent years there has been mounting evidence of the role played by magnetic fields 
at large scales [HI [T2] . Why should magnetic fields be assumed in various processes ranging 
from star formation to cluster dynamics and completely neglected prior to last scattering? 
Why are magnetic fields overlooked in CMB physics while they are observed in galaxies 
clusters, superclusters and high-redshift quasars? There are no reasons for doing so unless 
one would implicitly assume that large-scale magnetism suddenly arose between hydrogen 
recombination and, say, the gravitational collapse of the protogalaxy. While it might well 
be that the latter situation is the one preferred by nature, it would be nice to have some 
direct empirical evidence less biased by speculations. To comply with the latter program, a 
specific approach has been tailored through the last few years [T2] (see also [13]). The idea 
is, in a nutshell, to introduce consistently large-scale magnetic fields in all the steps leading 
to the estimate of CMB anisotropics and polarization. So far the program undertaken in 
[T2] led to various results 

• the large-scale magnetic fields have been included both at the level of the initial con- 
ditions as well as the level of the evolution equations for the standard adiabatic mode 
and for the other entropic initial conditions 



the temperature and polarization anisotropies induced by the magnetized (adiabatic 
and entropic) initial conditions have been computed 



the parameters of the magnetized background have been estimated (for the first time) 



2 A stands for the dark-energy component while CDM denotes the cold dark matter component. 
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by using the TT and TE correlationgj measured by the WMAP collaboration. 



There exist other approaches to the interplay between large-scale magnetic fields and CMB 
anisotropies (see [HI [T71 HH [191 120] for an incomplete list of references; see [12] for a more 
thorough account of earlier results). The common characteristic of those approaches has been 
to neglect the scalar modes of the geometry and to focus the attention to the tensor and 
vector modes. The recent results [131 [HI [15] show, in contrast with previous guesses, that 
large scale magnetic fields alter the initial conditions and the dynamics of the scalar modes 
of the geometry. They consequently distort, in a computable manner, the temperature and 
polarization anisotropies. 

A limitation common to nearly all studies on pre-decoupling magnetism has been so far 
the total absence of the effect of the magnetic field in the process of electron-photon scatter- 
ing. In [131 (HI US], for instance, the magnetic fields are included in the initial conditions and 
in all the relevant governing equations. The electron-photon scattering, however, is assumed 
to take place as if the magnetic fields were absent. The potential smallness of the effects does 
not justify its neglect since diverse small effects are often claimed to be detectable because 
of the purported control we now have on CMB foregrounds [21] . 

The consistent inclusion of magnetic fields in electron-photon scattering modifies quali- 
tatively the standard lore since the geodesies of electrons and ions are be modified by the 
presence of the Lorentz force term in curved backgrounds. Absent the contribution of the 
magnetic field, the motion of electrons and ions depends only upon the incident electric field; 
but when the magnetic field is included the classic treatment (see, for instance [I]) must be 
adapted to the new situation. 

The neglect of the role of the magnetic field in the electron-photon scattering has been 
recently relaxed in the guiding centre approximation [221 ESI [2U [22], and for a specified 
magnetic field orientation. The argument for keeping the direction fixed was essentially 
practical and in the present paper a general treatment will by developed along a twofold 
perspective 

• the orientation of the magnetic field will be kept arbitrary so that the matrix elements 
either in in the Jones or in the Mueller calculus will depend not only upon the directions 
of the incident and of the outgoing radiation but also on the magnetic field orientation; 

• after including the magnetic field in the scattering process the source terms for the 
transport equations of the scalar, vector and tensor modes of the geometry will be 
deduced explicitly. 

The latter analysis is still lacking both in the present and in the earlier literature. The 
magnetized electron-photon scattering is often required in diverse astrophyiscal situations 



3 Following the standard shorthand terminology the TT correlations denote the temperature autocorre- 
lations while the TE correlations denote the cross-correlation between the temperature and the E-mode 
polarization. 



like in the physics of magnetized sun spots [28J , or the theory of synchrotron emission 
whose results cannot be directly used since prior at last scattering electrons and ions are 
notoriously non-relativistic. Conversely some studies involving directly Thomson scattering 
in a magnetized environment [31] do not incorporate the fluctuations of the geometry and 
are also obtained using a preferential magnetic field orientation. 

The layout of this paper is therefore the following. In section [2] the tenets of the Mueller 
and Jones calculus will be reviewed and the matrix elements for the magnetized electron- 
photon scattering presented. Section [3] introduces the scalar, vector and tensor components 
of the brightness perturbations and the calculation of the collisionless part of the transport 
equations. The full scalar, vector and tensor transport equations will be discussed, respec- 
tively, in sections HJE3 and El Section [7] contains the concluding remarks. Explicit expressions 
involving all the relevant matrix elements both in the Jones and in the Mueller approaches 
have been collected in the appendices [A] and El 

2 Mueller and Jones calculus 

In the Mueller calculus the Stokes parameters are organized in a four- dimensional (Mueller) 
column vector whose components are exactly the four Stokes parameters, i.e. /, Q, U and 
V. In the Jones calculus the electric fields of the wave are organized in a two-dimensional 
column vector and the Stokes parameters are effectively derived quantities (see [27] for an 
introduction to the Mueller and Jones approaches). Hereunder a hybrid approach shall be 
employed. The polarization tensor "P%j = Vj% = E{ E* can be organized in a Stokes matrix 
whose explicit form is: 

= ^(I1 + Ua 1 + Va 2 + Qa 3 ), (2.1) 

where 1 denotes the identity matrix while ax, a 2 and a 3 are the three Pauli matrices. Some- 
times the Stokes matrix V is separated in a traceless part (i.e. the polarization matrix) 
supplemented by the identity matrix multiplying the intensity of the radiation field: this 
separation shall not be employed here. The orientation of the coordinate system is illus- 
trated in Fig. [TJ The radial, azimuthal and polar directions are 

f = (cos if sin ??, sin yj sin??, cos??), 

i? = (cos if cos??, sin if cost?, — sin??), 

tp — (— simp, cosy?, 0), (2.2) 

implying that fx'd — if. Photons propagate radially and n = (??, if) denotes the direction 
of the scattered photon while n' = (??', if') is the direction of the incoming photon; similarly 
H = cos?? and v = cost?'. 




When the photons impinge the electrons in a magnetized environment the magnetic field 
can be treated in the guiding centre approximation. Denoting with B the comoving magnetic 
field intensity the guiding centre approximation J32J [33] stipulates 

Bi(x, r) ~ Bi(x , r) + (x j - 4)9jBi + ... (2.3) 

where the ellipses stand for the higher orders in the gradients leading, both, to curvature and 
drift corrections which will be neglected in this investigation. The scales one must therefore 
compare are |x | = L , \x — Xq\ = L, A^ rec ' ) (the wavelength of the incident radiation at 
the recombination epoch) and H~l (i.e. the Hubble rate at recombination). It is easy to 
appreciate that A^ 60 -* = 0(/j,m) implying that 

#-* ~ L » L » Af c >. (2.4) 

Equation ( 12. 4 p implies that, for the purposes of the contribution of the magnetic field to 
the Stokes matrix the spatial gradients can be neglected while they cannot be neglected 
when estimating the effects of the large-scale inhomogeneities of the magnetic field. In 
spite of the fact that the contribution of the spatial gradients can be neglected in the first 
approximation, still the direction of the magnetic field should be appropriately taken into 
account. Consequently it is necessary to introduce a local basis which will define for us the 
magnetic field direction: 

ei = (cos a cos f3, sin a cos f3, — sin/3), 

e<i = (—sin a, cos a, 0), 

e 3 = (cos a sin (3, sin a sin (3, cos (3). (2.5) 

The basis of Eq. (12. 5p local since it accounts for the direction of the magnetic field over 
the typical scales involved in the electron-photon scattering. Once the direction of the local 
magnetic field has been fixed, the motion of the electrons and of the ions will follow the 
appropriate geodesies holding for charged particles in a gravitational field. 

For the calculation of the scattering matrix the magnetic field can be aligned along e%. 
The latter choice is purely conventional and it does not prevent from varying arbitrarily 
the direction of the magnetic field with respect either to the direction of propagation of 
the photons or to the direction of propagation of the other fluctuations of the geometry. 
Consider, as an example, a situation which will be treated later on in greater detail, i.e. the 
case where a relic vector mode of the geometry □ propagates along the direction k. Since k 
cooincides also, by definition, with the direction of the Fourier wavevector the whole problem 
will be characterized by 

• (n • k), i.e. the projection of the photon momentum along the direction of propagation 
of the relic vector; 



4 The same discussion, with the due differences, can be repeated in the case of the scalar or tensor modes 
of the geometry. Here the case of the vector modes is just selected for sake of illustration. 



• (e3 ■ k), i.e. the projection of the magnetic field direction along the direction of propa- 
gation of the relic vector. 

In the mentioned example we can choose, without loss of generality, k = z and the two 
physical polarizations of the relic vector will then be defined in the x — y plane. In this 
situation cos-$ = k ■ n and cos a = k ■ £3. The direction £3 does not coincide, in general, 
with z. For instance if a = (3 = — n/2, £3 coincides with e y while for a = and /3 = tt/2 
e% coincides with e x . This simple example shows explicitly that since the direction of £3 is 
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Figure 1: Schematic view of the relation between the coordinate system defining the scattered 
radiation field and the local frame of reference defining the direction of the magnetic field. 



arbitrary, the orientation of the magnetic field is also generic. Such an arbitrariness entails 
the dependence of the scattering matrix upon two supplementary angles. In total the Stokes 
matrix will then depend upon the two angles defining the direction of the scattered radiation, 
the two angles defining the direction of the incident radiation and the two angles defining the 
direction of the magnetic field. The Stokes matrix will then depend overall upon six angles: 
($, ip) (for the directions of the scattered photons), (i?', <p') (for the directions of the incident 
photons) and (a,/3) for the magnetic field direction. The schematic relation between the 
direction of the scattered radiation and the local frame defined by Eq. ( 12. 5ft is summarized 
in Fig. [TJ The (thick) dashed line denotes the direction of h, i.e. the direction of propagation 
of the radiation field. In Fig. [T]the two different alignments of e% are just meant to illustrate 
the effective arbitrariness of the magnetic field orientation. 

In the dipole approximation the scattered electric field can be computed as the compo- 
sition of the scattered electric fields due to the electrons and to the ions: 
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(2.6) 



where 0( e ) and am are, respectively, the accelerations for the electrons and for the ions. In 
the local frame defined by Eq. (12.51) the vector A = (o( e ) — <%)) the vector can be decomposed 
as A= (A^t + A 2 e 2 + A 3 e 3 ). Denoting with E\ = (E ■ ei), E 2 = (E ■ e 2 ) and E 3 = (E ■ e 2 ) 
the components of the electric fields of the incident radiation in the local frame we have that 
from the geodesies of electrons and ions 
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(2.7) 
(2.8) 
(2.9) 



where because of the global neutrality of the plasma, uq = fioa 3 is the common comoving 
concentration of electrons and ions; WBe,i and co pet i denote respectively the Larmor and 
plasma frequencies for electrons (and ions) 
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(2.10) 



where m e , ; denote either the electron or the ion mass depending upon the relative subscript 
and a(r) is the scale factor of a conformally flat geometry of Friedmann- Robertson- Walker 
type whose line element and metric tensor are defined as 



ds 2 = g^dx^dx" = a 2 (r)[dT 2 — dx 2 ], g^ v = a 2 (r)r) i 



fJIV ■ 



(2.11) 



In Eqs. (12. 7p . (12. 8p and (J2.9P the functions Aj (with i = 1, 2, 3) as well as ((lo) and all depend 
upon the angular frequency of the photon (i.e. u = 2-kv) and are defined as: 
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(2.12) 



The scale factor a(r) appears explicitly in Eqs. (I2.10p since the mass of the (non relativistic) 
species breaks the conformal invariance of the system of equations. Indeed, Eqs. ( 12. 7ft . 
(12. 8p and (12. 9p follow from the geodesies of charged species in the conformally flat metric 
of Eq. (12. lip where, for a generic massive particle, the mass shell condition implies that 
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g a pP a pP = m 2 (P a = mu a is the canonical momentum and u a the four- velocity) . Recalling 
that the comoving three-momentum q is defined as q = ap where Sijp 1 ^ = —gijP % P\ 
the comoving three- velocity is given by v = qj \/q 2 + m' 2 a 2 . Since the electrons are non- 
relativistic at last scattering q = mav and this is, ultimately, the rationale for the appearance 
of the scale factors in the explicit expressions of the Larmor and plasma frequencies for the 
electrons and for the ion^j. The numerical value of f e (co) for typical cosmological parameters 
is given by 

/.(") = (^) = 2.79 x 10- " (A) (2^) {Zm + 1)<<1} (213) 

where z* is the redshift to last scattering, i.e. z* = 1090.79lo!g2 according to the WMAP-7yr 
data [10J. In Eq. ( I2.13P B = \e 3 -B\; grossly speaking the typical values of v and B appearing 
in Eq. ( 12. 13ft do correspond, respectively, to the (very minimal) value of the frequency 
channel of CMB experiments and to the maximal value of the comoving magnetic field 
allowed by the distortions of the temperature autocorrelations and of the cross-correlations 
between temperature and polarization. The evolution of the Stokes matrix V can be formally 



written aqj 

dV 3f' r 

— + e'v = — m(q, q', a, p) v(n, n') M\n, n', a, /?), (2.14) 

dr 1d7t J 

where dQ' = d cos d' dip' and, defining the rate of electron-photon scattering T 7e , 

m, 



aT^ e = anox e a ei , a ie = -wr 2 , r e = — (2-15) 



is the differential optical depth. At the right of Eq. (j2.14p the matrix M(Q, Q', a, 0) is a 
2x2 and the four entries of the matrix M(Q, Q', a, 0) are separately reported in appendix lAl 
From the matrix elements of M(Q, Q' , a, (3) it is immediately possible to derive the evolution 
equations for the Stokes parameters in the Mueller form, namely, 

^ + e'X= J^ [dQ'T{Q,Q',a,f3)Z{Q'), (2.16) 

dr 62n J 

where X is a column matrix whose entries are, respectively, /, Q, U and V. The entries 
of the 4x4 Mueller matrix will be denoted as 71j(Q, Q', a, 0) where i and j run over the 
various Stokes parameters /, Q, U and V and are reported in appendix |A] (see, in particular, 
Eqs. (lA.lip - ( 1A.26P ). Finally, in terms of the matrix the evolution equations of the different 
Stokes parameters can be formally written as 



dl 

d~T + < 


s'l- 


dQ 
dr 


e'Q 



3e' 

32^ 
3e' 



32tt 



dQ'Fi(Q,Q',a,0), (2.17) 

JdQ'T Q (Q,Q',a,(3), (2.18) 



5 For more explicit discussions of these points see [TH1[S2]. Note that Eq. (A12) of [22^ contains few trivial 
typos which have been corrected in the archive version of the same paper. 

6 The dagger in Eq. (J2.14I) defines, as usual, the complex conjugate of the transposed matrix. 



ATT "if 1 i- 

-^ + e'U=—j dn'Fu(n,n',a,[3), (2.19) 

dV 3f' r 

-r + eV = — - / dn'jvcn, n', «, /?), (2.20) 

where in all the integrands at the right hand side of Eqs. (I2T7I) . (12~T8|) . (l2~T9j) and fl2T20|) the 
matrix elements 7ij(Q,Q',a, (3) are functions of, both, the angles of the incident radiation 
Q' = ($', f'), the angles of the scattered radiation Q = ("&,(p) and the orientation of the 
magnetic field defined by the angles a and (3. The explicit relations between %j(p., fi', a, (3) 
and the matrix elements Mij(Q, Q', a, 0) are reported in the appendix lAl 

3 Brightness perturbations 

The brightness perturbations, i.e. the fluctuations of the Stokes parameters in comparison 
to their equilibrium values can be decomposed as 

A x (x,r) = A? (x,t) + A£>(a?,r) + A?(x,r), (3.1) 

where X — I, Q, U, V denotes, generically, one of the four Stokes parameters and where 
the superscripts refer, respectively, to the scalar, vector and tensor modes of the geome- 
try. The scalar, vector and tensor components of the brightness perturbations are affected, 
respectively, by the scalar, vector and tensor inhomogeneties of the geometry and of the 
various sources. Assuming the conformally flat background introduced in Eq. (12. lip , the 
fluctuations of the metric can be written, in general terms, as 

5g^ u (x, t) = 5 s g^ u (x, r) + 5 v g^ u (x, r) + 5 t g^ u (x, r), (3.2) 

where S s , 5 V and 5 t denote the inhomogeneity preserving, separately, the scalar, vector and 
tensor nature of the fluctuations. The scalar modes of the geometry are parametrized in 
terms of four independent functions ipix, r), </>(x, r), E(x, r) and F(x, r): 

S s g 00 (x,T) = 2a 2 (r)0(f,r), 

8 s 9oi{x,r) = -a 2 (T)diF(x,T), 

S s 9ij(x, r) = 2a 2 (r)[^(f, r)(Jy - d.djEix, r)]. (3.3) 

By setting E and F to zero the gauge freedom is completely fixed and this choice pins down 
the longitudinal (or conformally Newtonian) gauge. The vector modes are described by two 
independent vectors Qi(x,r) and Wi(x,r) 



S v g 0i (x, t) = -a 2 Qi(x, r), S v gij(x, r) 



a 2 



d i W 1 (x ) T)+d 1 W i (x ) T) 



't"j \-^i 



(3.4) 



subjected to the conditions diQ % = and diW 1 = 0. It will be convenient, for the present 
purposes, to choose the gauge Qi = 0. The tensor modes of the geometry are parametrized 
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in terms of a rank-two tensor in three spatial dimensions, i.e. 

S t gij(x, t) = -a 2 hij, dih){x, r) = h\(x, r) = 0, (3.5) 

which is automatically invariant under infinitesimal coordinate transformations. The follow- 
ing shorthand notationL] will be adopted 

Cf\n, x, r) = d T A? + n%A? + e'A? + -(%), (3.6) 

q\dT/ s 

4 v )(n, x, r) = d T A { / ] + n%A^ + e'A^ + - ( ^L) , (3.7) 

q\dT/ v 

Cf\n,x,r) = d T A? +h%A? +e'A? + i(^) (3.8) 

where q = n^ 1 and where the scalar, vector and tensor contributions to the derivatives of 
the modulus of the comoving three-momentum are given, respectively, by 

^) = -qd T ij + qh%0, (3.9) 

p) =U, i n j (d i drW j + d T d j W i ), (3.10) 

arj v £ 

dq\ = _q 
dr) t 2 



^L) =-ln i n j d T h ij . (3.11) 



The identities of Eqs. (I3.9p . (I3.10p and (I3.1ip can be derived from the inhomogeneities of 
Eqs. (I3.3p -f l3~5|) by recalling the definition of comoving three momentum (see discussion 
after Eq. (12. lip ) and by using the relations 

where P % and P° are the space-like and time-like components of the canonical momentum 
obeying, for the photons, g a/ 3P a P 13 = 0. The notation introduced in Eqs. (13.61) . ( 13. 7ft . 
(13. 8p for the fluctuations of the intensity can also be generalized to the linear and circular 
polarizations: 

£ { P (n, x , r) = d T A { p + tidA^ + e'A { p , (3.13) 

where the subscript can coincide, alternatively, with Q, U and V (i.e. X = Q, U, V) and the 
superscript denotes the transformation properties of the given fluctuation (i.e. y = s, v, t). 
The fluctuations of the geometry would seem to affect only the brightness perturbation for 
the intensity but such a conclusion would be incorrect: in the presence of a magnetic field 
the evolution equations of the four brightness perturbations are all coupled by the collision 



7 The partial derivations with respect to r will be denotes by d T ; the partial derivations with respect to 
the spatial coordinates will be instead denoted by di with i = 1, 2, 3. 
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term which does not only contain the intensity of the radiation field but a weighted sum 
of the four brightness perturbations integrated over the directions of the incident radiation. 
Consequently the polarization of the metric fluctuations will also impact on all the four 
brightness perturbations. The conventions on the Fourier transform and polarizations of the 
scalar, vector and tensor modes will be, in short, 



(f)(x, r) 

W t (x,r) 
hij{x,r) 



(2tt) 3 / 2 
1 

0^/2 
1 

(2tt) 3 / 2 
1 



d 3 k (j)(k, t) e its , 
d 3 k^(k,r)e ils , 
d 3 kWi(k,r)e its , 

d 3 kh i:i {k,T)e ll£ , 



d i W\x,r) = 0, 
dihUx,T) = h\(x,r) 



(2vr) 3 / 2 
The vector and the tensor polarizations can be decomposed, respectively, as 

Wi(k,r) = J2e\ X) W ix) (k,r) = d t W a (k,T) + b l W b (k,r), 

A 

hij(k, r) = J2 4 A V) A r ) = e *%ft T ) + e®M£ t), 

A 



0. 



(3.14) 
(3.15) 
(3.16) 
(3.17) 

(3.18) 
(3.19) 



where k denotes the direction of propagation and the two orthogonal directions 6 and b 
are such that a x b = k. Supposing that the direction of propagation of the relic tensor is 
oriented along k, the two tensor polarizations are defined in terms of Oj and fej as: 



efJk) = didj - bibj, efJk) = dfij + dfii. 



(3.20) 



The projections of the vector and of the tensor polarizations on the direction of photon 
propagation n are: 



nW,(fc,r) 



h%W a {k,T) + n%W b (k,r) 



(3.21) 



Whijik, t) = U(h- a) 2 - (h ■ by]h @ (k, r) + 2(n -a)(n- b)h®(k, r) L (3.22) 

Choosing the direction of propagation of the relic vector and of the relic tensor along the z 
axis, the unit vectors a and b will coincide with the remaining two Cartesian directions and 
the related Fourier amplitudes will satisfy 



n¥i(fc,r) = 
h l h 3 hij{k, r) 




W L (k, t) Yf 1 ^ <p) - W R (k, t)Y?(&, <p) 
h R (k, t)Y 2 2 ($, <p) + h L (k, r)KT 2 (tf, <p) 



(3.23) 
(3.24) 
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where 



t W a (k,T)+iW b (k,r) t W a {k,r)-iW b {k,r) 
W L {k, t) = -= , W R (k, t) = 



h L (k,r) 



V2 
h®(k,r) +ih®(k,r) 



h R (k,r) 



V2 
h@(k,T) -ih®(k,r) 



V2 ' " V ' 7 y/2 

the spherical harmonics appearing in Eqs. ( 13. 23ft and ( 13.24ft are, respectively, 



(3.25) 



Y^0&,<p) = T\ — sin^ e ±Up 



¥?'{<&, ip) 



I15_ 
32tt 



sin 2 tf e ±2itp , 



(3.26) 



showing, as well known in the context of the total angular momentum method [34] , that the 
vector and tensor modes excite, respectively, the two harmonics given in Eq. (I3.26p . While 
the total angular momentum method can be generalized to the case of an arbitrarily oriented 
magnetic field, we prefer to work, in the present context, with the formalism which is more 
directly applicable to numerical codes and to standard analytic estimates. In Fourier space 
Eqs. (13TB]) . (1577) and (EH) become 



C { j\fi,if,k,r) 
Cy(ji,(p,k,r) 

cf\fj,,(p,k,T) 



+ 



d T Af } + (ik/i + e') Af + ik^cj) - d T ip, 

drA^ + iikfi + ^A^ 




I fj, 



d T W L (k, r) Yf 1 ^, <f) - d T W R (k, t)Y?{#, <p) 




!t) + {ik^i + e')Af 
d T h R (k, t)Y 2 2 ($, tp) + d T h L (k, r)Kf 2 (tf, tp) 



Similarly Eq. ( 13.131) becomes, in Fourier space, 

4 y V tp, k, r) = d T A { p + {Hep + e')A%\ 



(3.27) 
(3.28) 

(3.29) 
(3.30) 



The explicit form of the transport equations for the scalar, vector and tensor modes of the 
geometry will be scrutinized in the three forthcoming sections. 



4 Scalar modes 



Following the notation of Eqs. f 13 . 2 Tj) and ( I3.30p the scalar transport equations can be 
formally expressed as 



<-) - ( i 3e' f 1 f 2lT ( ) 

£/ ; (/z, tp, k, r) = e'h l vf> + — - du dip'J 7 } 8 ' (/i, u, <p, tp', a, 0), 

6lix J-\ Jo 

<-) - 3e' r 1 r 2n i ) 
£q (/wA^) = — J * du J (i^ !, (/i,y,^^a,|3), 



(4.1) 
(4.2) 
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Cff (//, if, k, t) = — - / dv / V J# ; (//, i/, </?, y/, a, /3) , (4.3) 

OZ7T 7-1 JO 

Cv{v,V,k,T) = —— dv I d(p' J^ } (ji,v,<p,<p',a,fi), (4.4) 

OZ7T J-l JO 

where u $ - denotes the scalar component of the baryon velocity field. In Eqs. ( I4.ip - (j4.4p 
the source terms involve the integration over the incoming photon directions. Both the 
integration over ip' and v can be performed explicitly and the final expressions are rather 
lengthy, as easily imaginable. To make the explicit equations more manageable without 
loosing any relevant information it is useful, in the following part of the present section, 
to write the results already in the physical limit, i.e. owing to the numerical values of the 
plasma and Larmor frequencies and recalling Eq. (12.121) 

Ai(w) = A 2 (w) = A 3 (w) = 1 + 0{m c /m p ), 

CH^-l + / e 2 H + £>(/"*). (4.5) 

The scalar source terms depend upon the explicit form of the matrix elements appearing in 
Eqs. (!A.7p . (lA.8p . (1A.9J) and (lA.lOp . The integration over ip' can be performed explicitly. 
Using the notation 

Tab(v, v, V, a, p) = / dtp' Tahiti, v, (p, <f/, a,P), (4.6) 

JO 

the final results are reported, for completeness and future peruse, in appendix [Bj According 
to Eqs. (I4.ip - P~4"p . the expressions reported in Eqs. ( ]B.ip -( )B~T6|) must be integrated over v. 



For the v integration it is useful to expand the various brightness perturbations in a series 
of Legendre polynomials Ptiy) 

A x (u, k, r) = Y,(-i) l W + 1) Pl(") A xe(k, r). (4.7) 

i 

The integration over v will then have the net result of expressing the source terms in terms of 
a limited number of multipoles of the intensity and of the polarization. In explicit terms the 
source terms can be expressed, for each brightness perturbation, as an expansion in f e {w): 

d T A® + (ik/i + e') Af = d T ^j - iky4 + e'Ai + d /,(w) B T + d f c {u) C h (4.8) 

d T A$ + (ikti + d) A§ } = dAq + d / a (w) B Q + d /» C Q , (4.9) 

d T A ( fi + {ikft + e)Aff = e'Au + d / e (w) B v + d ft (w) C v , (4.10) 

d T A$ + (ikti + d)A$ = e'Av + d f c (u) B v + d f 2 c {u) C v , (4.11) 

where, for X = I, Q, U, V, Ax denotes the leading order result, Bx denotes the next- 
to-leading order (NLO) correction while Cx denotes the next-to-next-to-leading (NNLO) 
term. Defining with Sp the usual combination of the quadrupole of the intensity and of the 
monopole and quadrupole of the linear polarization (see, e.g. [21 El H]) 

S P = A I2 + Aq + A Q2 , (4.12) 
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the leading order contribution for the for brightness perturbations is then given by: 

Ai = A I0 + fiv h -^-S P , (4.13) 

Aq = 1(1-H 2 )Sp, (4.14) 

3 

Au = 0, A v = --ifiAvu (4-15) 

where the notation ^ = kv\, has been employed for the scalar component of the Doppler 
term. The NLO contribution to the right hand side of Eqs. (|Ojl . ( 1431) . (OPjl and (1431"]) is 

3 



Bq 
Bu 

By 



iA 



1 + /i 2 ) cos j3 + fj,J 1 — /i 2 cos (<y9 — a) sin /3 A 
va (/i 2 — 1) cos/3 + /W 1 — /i 2 cos (<p — a) sin/3 



VI 



- i Ayi y 1 — /i 2 sin /3 sin (<p — a) , 

Vi -/i 2 . ■ 

/i cos p cos (<£> — a) sin p 



/l 2" 

2/i cos (cp — a) sin 2 /3 H sin /3 



(4.16) 
(4.17) 
(4.18) 



— cos p cos [if — a) sin p 



A 70 

A/2 

(A Q2 + Aq ) 



(4.19) 



Finally, the NNLO contribution to the right hand side of Eqs. (jMD, (Oj) . ( H1D]I and (gXQ) 

is 

•A 70 



C; 



(// + 1) + /iv/l — /i 2 cos (cp — a) sin 2/3 



A 
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1 

+ 2 



ii 2 (A /0 + A J2 ) + (A 72 - 2A /0 ; 

l + /x 2r 



n 9 / \ CUb "/-' / 9\ 

sin pcos {ip — a) — (1 + 3/i ) 



|(1 - fi 2 ) sin 2 (p - a) + 
\i\J\ - /i 2 



cos 2(^9 — a — /3) + cos 2(<^ — a + /3) 



sin (<£> — a — 2/3) — sin (cp — a + 2/3) 



C, 



Q 



/i 2 - 1 sin 2 /3 



/i 2 - 1 sin 2 /3 

1 - 

4 2 

,2 1 ,,2 



4 - 2(2/j 2 + 1) cos 2 (<p - a) 



+ 



(Aq + Aq 2 ), 



(4.20) 
cos (<^ — a) sin 2/3 > A/ 



1 + (1 — /i 2 ) cos 2 (if — a) 



H cos (<£ — a) sin 2/3 }> A/ 2 



i^—l-^-t±cos2((p-a) + l(l-fi 2 )cos2P 



+ 
+ 



/i 2 + l 
16 

/Vl - A* 2 
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cos 2(ip — a — (3) + cos 2((p — a + (3) 
sin {ip — a — 2/3) — sin (99 — a + 2/3) 
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(A Q0 + A Q2 ), 



(4.21) 



sin {ip — a) 



VT^JF 



sin 2/3 — fi cos (ip — a) sin 2 (3 



+ sin {ip — a 

sin (3 sin {ip — a 



fi cos {ip — a) sin 2 f3 sin 2/3 



A/2 



+ 



Cv 



3 

--2 

2 



2 
/j- cos 2 /3 



1 — ji 2 cos /3 + 2/i cos (<£> — a) sin /3 

Ayi. 



(A Q2 + A Q0 ), 



/i 2 cos (</? — a) sin /3 cos /3 



(4.22) 
(4.23) 



Several cross-checks on the obtained results have been made; they will be swiftly mentioned 
and can be directly reproduced by using the results reported in the appendices [A] and El 

• it has been verified explicitly at the level of the exact expressions (i.e. without imple- 
menting the limit of Eq. ( |4.5p ) the equations must be independent upon the a and /3 
once (£3 • B) — > 0: this is exactly what happens; 

• it has been verified that in the limit a = (3 = the exact expressions must reproduce 
the partial results already obtained in [231 I2S1 [26] ; with Eqs. ( 14 . 8 [) ( 147TTJ) few typos 
present in the published version of [23] are corrected; 



• by averaging of the source terms over a and /3 terms proportional to f e {oj) should 
automatically disappear without performing any specific limit: this is what will be 
explicitly shown in the remaining part of this section. 

The remaining part of the section is devoted to the averaging of the source terms over 
the magnetic field directions as suggested in the last point of the above list of items. By 
integrating over a and (3 the source functions appearing at the right hand side of Eqs. ( 14. ip . 
( 14. 2p . ( 14. 3 p and (14. 4p . the evolution equations for the brightness perturbations read 



.(8), 



e'nV (s) 



Cy(ji,<p,k,T) 

3c' /"i r^ n r 71 r^ n 1 \ 

+ _^ „ / dv / dp / sin/3<i/3 / daT] (/1, z/, ip,ip', a,/3), 
J-i Jo ' Jo Jo 

' [ • ' '" dip' [" sin /3d/3 f " daJF§{pi, v, <p, tp\ a, /3),(4.25) 



1287T 2 7-1 

^n(/i,^,fc,r) 



(4.24) 



-Q 
.(b) 



£\j'{fi,ip,k,T) 

C { v } {fi,<p,k,T) 



, dv 

1287T 2 7-1 Ju 

3e' 



128tt 2 7-1 
3e' 



2 k 



:W 



sm j3d/3 / daFv {fi,i>,tp,tp', a, (3), {4.2G) 



1287T 2 



1 i-2-k 

dv \ dip' 

Jo ' Jo Jo 

dv f ^ dip' [" sin. Pdfi I ^ da^ s> {fi, v, ip, <p', a, /3), (4.27) 
1 Jo Jo Jo 



A*), 



where the factor 1287T 2 accounts for the 47r factor arising from the average over the solid 
angle spanned by a and (3. By performing the averages explicitly, the evolution equations 
of the four brightness perturbations read: 



» 



d T Af> + {ikfi + e')A ( i 



f\ A (s) 



d T ip — ikfi(j) + e' 
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A/o + [iv h -— S P 





+ <|a. + ^*)] 


d T A§ + (ikfi + e') A§> 


,(/e 2 -3)(/i 2 -l), 

= e 4 6 P , 


d T A® + (ikn + (T)A$ 


= o, 


d T AP + {ikv + e') A$ 


= -^(3 + / e 2 )A vl . 



(4.28) 

(4.29) 
(4.30) 
(4.31) 

The results of the present section make quantitatively clear that if the magnetic field has a 
predominant direction over typical scales comparable with the wavelengths of the scattered 
photons, then the circular polarization is larger than in the case where, over the same physical 
scales the magnetic field is randomly oriented. It has been argued in [22], [23] (see also [26] ) 
that over small angular scales the maximal amount of circular polarization arises when there 
is a strong alignment of the magnetic field along the direction of propagation of the photon 
which coincides, for large multipoles, with the third Cartesian direction. 

The present results improve and confirm, at once, the assumptions made in the analytic 
and numerical estimates of magnetized CMB anisotropics of Refs. [T3l [T4"| [15]. Indeed, 
the scalar fluctuations of the geometry obey a set of evolution equations where large-scale 
magnetic fields contribute in many respects. These equations will not be repeated here and 
can be found in p3] [HI US]- The present results improve on the transport equations used 
there and pave the way for a more consistent account the effects of pre- decoupling magnetic 
fields both at large as well as at small angular scales. 

5 Vector modes 

In the case of the vector modes of the geometry the integration over <p' of the source func- 
tions cannot be easily performed as in the case of the scalar modes of the geometry (see 
section H]). Each of the two vector polarizations induce a different angular dependence in the 
corresponding brightness perturbations. In this paper three categories of circular and linear 
polarizations can be defined: 

• the linear and circular polarizations of the scattered (and incident) photons (already 
described in sections [5] and [3]) which are described by the four Stokes parameters or by 
the appropriate Stokes matrix; 

• the linear and circular polarizations of the relic vector waves which we are going to 
discuss in the present section and which have been already introduced, respectively, in 
Eqs. fl3~2T]) and (13T23J) ; 

• the linear and circular polarizations of the relic tensor waves introduced, respectively, 
in Eqs. fj3.22j) and (I3.24[) and discussed in the following section [U 
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The linear and circular polarizations of the relic tensor and vector waves are just equivalent 
basis for the description of the tensor and vector modes of the geometry. To avoid potential 
confusions the vector and the tensor waves will always be treated in the basis of the linear 
polarizations. In full analogy with the treatment of section H] the evolution equations for the 
vector components of the brightness perturbations can be formally written as 

4 v W,£,r) = e ' n \ (v) + £L f 1 d v r dp'^ v \n,v,p,p',a,P), (5.1) 

OZ7T J-l JO 

( ) - 3 e ' f 1 f 27r ( ) 

£q (fi,<P,Kr) = ^-y_ d » ) dp' Fq (ji,v,<P,<p',<x,P), (5.2) 

'( v )^,, ,„ Z ^\ - 6e / A,. / A, J X-(v)/ 



C\J'{n,p,k,r) = — j dv J dtp' J$ } (fi,u,(p,(p',a,p), (5.3) 

£<v(p,'P,k,T) = — dv I dtp' jty } (ii,i>,<p,<p',a,0), (5.4) 

OZ7T J-l JO 

where v | v denotes the vector component of the baryon velocity field. The polarizations of 
the baryon velocity will follow the same kind of decomposition illustrated for the vector of 
the geometry in Eq. (13. 23ft . The relative directions of the magnetic field intensity and of 
the photon propagation determine the polarization of the outgoing radiation. Following the 
strategy described in section [2] the direction of propagation of the relic vector wave can be 
fixed and the direction of the magnetic field varied at wish. 

Consider first the case where the magnetic field is oriented along the same direction of 
the vector wave and suppose, without loss of generality, that the vector propagates along 
k = z. Since, in this case, a = (3 = and Eqs. (lA.ljl . (IA.2I) . (1A.3I) and (IA.4I) will lead, 
respectively, to the following matrix elements: 



M u (fi,p,v,p') = C/iz/Ax cos (<p' - p) — \] 1 - ^a/I - v 2 A 3 

- iA 2 fcCfivsm (ip'-p), 

Mi 2 {n,p,v,p') = -C/iAisinA^ - iA 2 f e ([icos(p' -p), 

M 2 i(n,p } v } p') = (vA ± sin Ap + if c A 2 (v cos (p' - p), 

M 22 (n,p>,v,p>') = CAi cos {iff- ip) - if C A 2 C, sin (p>' - p). (5.5) 

Using Eq. (15. 5ft the source terms of Eqs. ( 15.11) . (15.21) . (15. 3ft and (15.41) can be computed. 
The explicit form of Eqs. (I5.1I) - (I5.4I) is rather lengthy: instead of writing all the equations, 
Eq. (15.11) will just be written for illustration with the purpose of demonstrating how the 
different vector polarizations induce a specific azimuthal dependence in the vector brightness 
perturbations. Equation ( 15.11) written in the basis of the linear polarizations reads 



W _l (n.,, _l c 'uW _i_„- „ . /1 _ „2 



d T Ay> + (ikfi + e')Ay> +i^Jl-^ 



cos <pd T W a + sin <pd T W\ 



3e' f l , r 2 ^ 



\j\ -/i 2 [cos^ a + sinv9t; fe ] + — - / dv dp)' J} vj (/x, v, <p, ip', a,fl), (5.6) 



327T J-i Jo 
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where, according to Eq. ( \A.7\i the integrand of the source term acts on the vector components 
of the various brightness perturbations and it is given by 

jj v \fi, is, if, <p', a, (3) = Tii(fi, is, (p, (p', a, f3)A { /\v } ip') + Ti Q (fi, u, <p, tp', a, /3)A ( q\u, ip') 



+Tiu(p, ", <P, <p', ®, P) A u\v, p') + 7/y(/i, v, (p, <p', a, l3)A ( y\is, ¥>')■ 



(5.7) 



After inspection of all the four expressions appearing in Eqs. (15. ip . (15. 2p . (15. 3p and (15. 4p it 
can be checked that the consistent ansatz for the four brightness perturbations is given by 

A^ ] (<p,n,k,T) = 



1-f 



&V \(p,/JL,k,T 

A ( y\ip,fi,T) = 



cos ipM. a {k, r) + sin (pAib(k, r) 

cos ipj\f a (k, r) + sin cpHb{k, r) 
• sin pj\f a (k,r) + cospAfb(k, r 



\-p? 

cos pV a (k, t) + sin pVb{k, r 



1 — fi 



(5.8) 

(5.9) 

(5.10) 

(5.11) 



The equations obeyed by A4 a , N a and V a are the same as the ones obeyed by Ai b , J\f b and 
Vb and they can be written, for a generic linear polarization, as 

d T M + (ikfi + e')M + ipid T W = e'v + eVCAiA 3 s[ v) + e / /i/cCA 2 A 3 S^ ) , 



d T M + (ifc/x + e')M = e'CAiA 3 sS v) + e7 c A 2 A 3 E< v) , 



d T V + (ikfi + e')V = e7eCA 2 A 3 X (v) + e'CA^E^, 

where the two newly defined source functions E^ and E 2 are given by: 

3 n 



(5.12) 
(5.13) 
(5.14) 



i-i 

3 rl 



,u J u 1 „: ,YI i ^ (v) 

.2 i\iiCL ,, _\ 



is(is 2 - l)M(k, u, r) + (z/ 4 - l)Af(k, is, r) 



dv, 



(5.15) 
(5.16) 



^ ( i\k,r) 

X%\k,T) = ^J y-l)V(k,u,r)du. 

The source functions appearing in Eqs. (15.15P and (I5.16P can be made more explicit by 
expanding M.{is, k, r), M{is, k, r) and V(is, k, r) with the same conventions employed in Eq. 
@2J: " 

M(is,k,r) = ^2(-i) i (2£+l)P l (u)Mi(k,r). 

e 

M{u,k,r) = £H)^+1)P^)A/KM- 

e. 

V{v,k,r) = £H)'(2*+1)P/0>)V/(A;,t), 



(5.17) 
(5.18) 
(5.19) 



where, following the same conventions of Eq. (14. 7p . M.^, Ng and Vi denote the £-th multipole 
of the corresponding quantity. The result of the integration over v is therefore 

3 



sfW) = ^-|*i-^ 
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i(Mi +Mi 



SV'lt.r) = -£-£, 



^t). 



(5.20) 
(5.21) 



The same considerations developed in the basis of the linear vector polarizations can be 
repeated in the case of the left and right polarized waves. Bearing in mind Eq. (I3.26p . Eqs. 
(I5T8]) -( I5TL|) can be written: 



VI 



A\P((p,fj,,k,r) = -2% 



2// 



Yf 1 ^, <p)M L (k, r) - Yr\ti, <f)M R (k, r) 
Yf 1 (n,<f,k, r)Af L (k, r) - Yf 1 (ji, <p)Af R (k, r) 
Yf 1 ^, <f)Af L (k, r) + Yr\», <f, k, r)Af R (k, r) 



A{J } (<p, im, k, t) = iT- Y-\n, v )V L (k, t) - YfiQi, v)V R (k, t) 



(5.22) 
(5.23) 
(5.24) 
(5.25) 



For sufficiently small angular scales (i.e. for sufficiently large multipoles) the microwave sky 
degenerates into a plane and the. In this situation microwave photons propagate, for all 
practical purposes, along the z axis and instead of the spherical decomposition based on 
spherical harmonics one can safely use a plane-wave decomposition. Since the wavelength of 
the photons is typically much shorter than the inhomogeneity scale of the magnetic field one 
could also argue, at this point, that the situation in which the magnetic field oriented along 
the direction of propagation of the relic (vector) wave is sufficiently generic. This conclu- 
sion should however be scrutinized more carefully and this is the purpose of the discussion 
reported hereunder. 

Suppose that the direction of propagation of the vector wave is not parallel to the mag- 
netic field direction but orthogonal. If the relic vector propagates along the magnetic field 
direction, then, in the language of Eq. (12. 5p . k || £3 implying a = (3 = 0. If the relic vector 
propagates orthogonally to the magnetic field direction then we can set a = (3 = —n/2 
implying that k _L £3. The direction of k will still be chosen to be the z axis so that 
(k ■ h) = n = costf. In the case a = (3 = -n/2 Eqs. flA.lj) . (IA.2J) . (IA.3J) and ( 1A.4J) read 



M n (fJ.,<p, v,(p') 
M 12 (/j,,(p,v,(p') 

M 22 (n,(p,v,ip') 



(Aji-^VT 



CAifj,u cos ip' cos ip — A 3 /iz/ sin ip sin ip' 



+ if c A 2 ((i , yl — /i 2 cos p>' — fiyl — v 2 cos <p) 
= — (fiAi cos ip simp' — A 3 fi cos (p' simp 
- if e (A 2 Jl- jj? sin ip', 



—A 3 u cos ip sin cp' — CAiZ/ cos (p' sin ip + if e (A 2 Vl — v 2 -, 
£Ai sin ip' sin ip — A3 cos ip' cos ip. 



(5.26) 

(5.27) 
(5.28) 
(5.29) 



In this case we can already expect, in comparison with the situation k \\ e%, that the transport 
equations differ depending upon the specific vector polarization. The solution of the system 
can indeed be written as 



.(v) 



Ay(ip,fi,k,r) 



p 



cos ipM. a (k, t) + sin ipMb(k, r) 



(5.30) 
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^Q (<P, V, k, T) = lly/l - /I 2 

A { y\<p,n,k,T) = J\ -/i 2 sin2^V a (/c,r) + fiV b (k,r 



l-/i 2 



cos ipj\f a (k, t) + sin ifj\f b (k, r) 
■ sin <fj\f a (k, r) + cos ipAf b (k, r 



(5.31) 
(5.32) 
(5.33) 



The azimuthal factorization of Eqs. (I5.30p - fl5.33p is not arbitrary and it is dictated by the 
specific form of the system of Eqs. (I5.ip - fl5.4p in the case when k _L £3. Starting with the 
polarization W a the corresponding evolution equations for Ai a , M a and V a are given by: 



d T M a + {ikn + e')M a + i»d T W a - e'v a = -e'v( 2 (AJ - f 2 A 2 2 )^, 
d T Ma + (ikiM + e')M a = -e'K 2 (Ai - /eA^sS, 
d T V a + (ik/i + e')V a = 0, 



where 



V (v) 



du[u(u 2 -l)M a + (u 4 -l)K}. 



(5.34) 
(5.35) 
(5.36) 

(5.37) 



Consider then the case of the polarization W b . The evolution equations are, in this second 



case, 



d T M b + (ikfj, + e')M b + ifj,d T W b - e'v b = e'/iCAiA 3 £ 



(v) 
16 



(v) 
16 



-~/ e CA 2 A 3) u j\ 2 V h 

d T Af b + {ikn + e')Af b = e'CAiAsS 

~f e (A 2 A 3 j\ 2 V b dv, 

d T V b + (ikfi + e)V b = e7 e /iCA 2 A 3 S 

3 f 1 

--/eCAxAg/iJ v 2 V b dv, 



(v) 
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where 



.(v) 

J 16 



-1 



dls[ls(ls 2 -l)M b +(ls 4 -l)N b }. 



(5.38) 



(5.39) 



(5.40) 



(5.41) 



The two sets of equations reported in Eqs. f l5.34p - fl5.36l) and in Eqs. f l5.37p - fl5.40p show 
various interesting features which can be summarized as follows: 

• if e 3 || k (i.e. a = (3 = 0) the evolution equations of the two vector polarizations are 
independent insofar as they can be given different initial conditions prior to decoupling 
but the evolution equations of the corresponding brightness perturbations are the same; 

• if £3 _L k (e.g. a — (3 — — tt/2) the two linear polarizations are equally independent 
but obey different evolution equations as it is clear by comparing Eqs. (15.340 - 05. 36D 
with Eqs. fl53B]) - fl5^0]) : 
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• the equation for the a- vector polarization (i.e. Eqs. f )5.34p -f l5.36p ) lead to linear and 
circular photon polarizations which are a factor 0(f e ) smaller than the corresponding 
equations for the b- vector polarization (see Eqs. f)5.37p - fl5.40p ). 

When k _L £3 and k = z we also have that a = x and b = y. But if a = fi — 
— 7r/2, then £3 = y. Therefore the amount of magnetically induced linear and circular 
photon polarization is larger when the magnetic field and the vector polarization are oriented 
along the same direction. In Fig. [2] where the geometric set-up of the vector problem 
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Figure 2: The interplay between the two linear vector polarizations (i.e. a and b), the 
magnetic field direction and the direction of propagation of the scattered radiation (i.e. h). 

is summarized. The wiggly line represents pictorially a vector wave propagating in the 
direction k which has been taken to be aligned with the z axis. Always in Fig. [2] the shaded 
plane denotes the polarization plane of the vector wave spanned by the two unit vectors a 
and b. Finally h denotes the direction of propagation of the photons. If the direction of the 
magnetic field is parallel to the direction in which the vector modes propagate (thick arrow 
in Fig. [2]), the photons do not inherit a computable amount of circular polarization and, 
furthermore, the two linear vector polarizations will lead to the same transport equations 
for the brightness perturbations. Conversely, if the magnetic field is parallel to one of the 
two vector polarizations (thick dashed arrows in Fig. [2]) the transport equations for the two 
linear vector polarizations will be different. If the linear vector polarization is aligned with 
the magnetic field intensity (for instance B || b) the transport equations for Ai a , M a and 
V a will lead to a V-mode polarization larger than the one generated by the other vector 
polarization and described in terms of A4&, A4 and V&. 

The source terms for the evolution equations of the vector modes can be averaged over 
the orientations of the magnetic field, i.e. 

£/ v) (/i,y?, k,r) 



e'tivV 
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1287T 2 J- 



1 r2ir rn r2-K , . 

dv \ dip' / sin (3d(3 \ daJ-'j (jj,,v,ip,ip',a, /?), 
-1 jo jo jo 



(5.42) 



, . _, gg' /-l /■27T /-7T rlis . . 

£q (v,¥,k,T) = — — I dv I dip' / sin /3d/3 / da J 7 ^ (fi,v,<p,<p',a, 13) ,(5.43) 

12071"^ J-l JO JO JO 

, . _, Qzr' ("1 /■27T /"JT /■27T , . 



Hj (v,<P,k, T ) 



,(v) 



£y } (fi,(p,k,T) 



128tt 2 y_i 
3e' ^ 



1287T 2 



/•27T /*7T /•27T 

di/ / dcp' / sin/3<i/3 / da.Fy^/i, z/, </>,<//, a, /3). (5.45) 

1 JO JO JO 



r(v), 



The direct computation of the averaged source terms leads to the same expression for both 
vector polarizations. Denoting with M. either Ai a or hA h (and similarly for M and V), Eqs. 
( I5.42p . (I5.43p . (15 .440 and (I5.45J) lead to the following triplet of equations: 



d T M + (ikfi + e')M + iftd T W - e'v = —e' 

15 



$m - 7A? 



+ 6CAiA 3 -2A^ 



n 



(v) 



d T M + {ikfi + e')M 



d T V + {ikfi + e')V 



t.4 

15 



5/ ^ - 7A^ 



+ 6CA!A 3 -2A? 



Si 



(v) 



(5.46) 



(5.47) 



r C[2A!A3-C(Ai + /^)]^. (5.48 



,(v) 



As in the case of Eqs. (I4.28p - fl4.3ip if the magnetic field has a predominant direction 
over typical scales comparable with the wavelengths of the scattered photons, the circular 
polarization of the photons induced by the vector modes is larger than in the case where, 
over the same physical scales the magnetic field does not have a specific orientation. This 
conclusion can be reached by comparing Eqs. ( 15.46(1 . (15.471) and (" 15.48(1 to Eqs. (I5.37p -f l5.40p 
obtained in the case when the magnetic field is oriented along one of the two polarizations 
of the relic vector. 



6 Tensor modes 



Recalling the notations introduced in Eqs. fl3.27p -f l3.30p . the evolution equations for the 
tensor components of the brightness perturbations shall be written, in general terms, as 



'Xf 1 rl /-27T 



Cy(ji,(p,k,r) 



32tt J -i 
3e' r 1 



d<P , J r i t) (lJ>,v,<P,<p',(x,P), 



32tt j"-i 
3e' r 1 



2n 



(li > I dip' J^'(ji,v,(p,(p',a,^), 



2n 



<*), 



dv dtp' J^'(jj,,v,(p,(p',a,P), 

62,71 J-i Jo 

3e' r 1 r 27T (+\ 

— — / dv dyj' jy(/x, z/,y?,y?',a,/3). 

3Z7T J-l JO 



(6.1) 
(6.2) 
(6.3) 
(6.4) 
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Consider first the case where the propagation of the relic graviton is parallel to the direction 
of the magnetic field intensity. The direction of propagation of the tensor wave can be 
chosen, without loss of generality, as k = z. Therefore Eq. ( 12. 5 j) implies that a = (3 = 0, i.e. 
k || £3. To illustrate the azimuthal dependence of the problem it is instructive to write down 
Eq. (16. ip in explicit terms: 



3e' /■' 



.(*) 



1 



d T Ay > + {ikn + e')A?'--{l -,r 



cos 2ipd T h S) + sin 2ipd T h i 



dv \ dip' ff ] ((j,,v,(p,(p',a,P), 



r ., , , r . . . r ,. , (6-5) 

62n J-i Jo 

where, according to Eq. (IA.7J) the integrand of the source term acts on the vector components 
of the various brightness perturbations and it is given by 

■7/ (p, ", <P, <P', a, P) = 7//(/x, v, <p, if', a, /?) Af } (v, ip') + Tiq(ji, v, <p, <p', a, fi)A% (u, ip') 
+Tiu{n, ", <P, v', ®, /3)A|) ) ( z/ > f') + Tivil-i, v, <p, (p', a, (3)A ( y\v, ip'). (6.6) 



The remaining three equations (i.e. Eqs. (16. 2p . (16. 3p and (16. 4p ) have a similar structure but 
the contribution of the tensor modes of the geometry is absent. The azimuthal dependence 
can be decoupled from the radial dependence and the brightness perturbations will be 



A ( j\ip,fi,k,T) 
A^(^,/i,/c,r) 
A§\(p,(i,h,T) 
A$((p,(i,h,T) 



fl") cos2ipZ e (/j, } k, t) + sin2y2,2^(/i, k, r 
(1 +/i 2 ) cos2ipT (S ([i,k,T) + sin 2ipT® (n, k, r) 
2/i — sin 2ipT (S) {[i, k, r) + cos 2ipT®{ix, k, r 
2/i cos 2ipS (B (/i, k, r) + sin 2ipS® (/i, k, r) 



(6.7) 

(6.8) 

(6.9) 

(6.10) 



In the case k \\ e% the symmetry of the system implies necessarily an ansatz in the form of Eqs. 
(16. 7p . (16.81) . (16.91) and (I6.10p . Starting with the explicit form of Eq. (16. 5p . it is immediately 
clear that the form of A\' (y?, /jl, k,r) is constrained by the ip dependence appearing at the 
left hand side of Eq. (16. 5p . The integrand at the right hand side of Eq. (I6.5p . i.e. Eq. (16. 6p . 
contains also Aq (ip 1 , u, k, r) whose explicit form is univocally determined by observing that 
the integral over ip' must match with the ip dependence appearing at the left hand side of Eq. 
(16. 5p . But the obtained ansatz for Aj(ip,fi,k,r) and Aq ((p,fi, k, r) can be inserted back 
into Eq. (16. 2p : this step will constructively determine the explicit form of A n (y3,/i, k, r). 
Equation (16. 3p will finally determine the explicit form of Ay (ip, /i, fc, r) whose ip dependence 
will have to be consistent with Eq. ( 16 .4p . The result of this procedure, expressed by Eqs. 
(16. 7p . (16. 8p . (I6.9P and (I6.10p . determines the radial evolution and, in particular, the following 
set of equations [25] : 



d T Z + (ikfx + e')Z — -d T h 



e'( 2 (L;)[Al(u) - fMAKu)}^, 



(6.11) 
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d T T+(zkfi + e')T+e'T=-e , C 2 (uj)[Al(uj)-f^uj)Al(uj)}^ t \ (6.12) 

d T S + (ikti + (?)S = Q, (6.13) 

where Z, T and S denote either the © or the © polarization. By expanding Z, T in series 
of Legendre polynomials 



Z(u, k 7 t) = £H)'(2£ + 1) P t (u) Z t (k, r). 
l 

7>, k, r) = J2(-iY( 2i + 1) ^0) TK^ r), 



(6.14) 
(6.15) 



the source term E^ can also be expressed as 

3 
32 7 

-2^2 -2o 



S (t) = ^ £*/[(!- i^) 2 Z(i/)-(1 + 1^)^(1/) -4i^T(i/)] 



70 Z4 + T"To"70 r4+ 7 r2_ 5 To ' 



(6.16) 



where, as usual, Zi and 7^ denote the £-th mulipoles of the corresponding functions. The 
results obtained in Eqs. (16. lip . f)6.12p and (I6.13P with the partial treatment of the tensor 
modes developed in [25] • As in the case of the vectors instead of working with the linear 




Figure 3: The © and © polarizations are illustrated, respectively, with full and dashed lines. 
The direction of propagation of the wave is not shown and it is orthogonal to the plane 
spanned by a and b. When B || k the magnetic field is oriented perpendicularly to the plane 
of the figure. 

polarizations of the relic gravitons, we could as well work with the circular polarization. 
The results obtained so far can be easily translated to the case when the relic gravitons are 
circularly polarized, always assuming that k \\ e%, i.e. that the direction of propagation of 
the relic gravitons is parallel to the orientation of the magnetic field intensity. 
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As expected from the vector case, when the relic tensor propagates orthogonally to the 
magnetic field direction, the two tensor polarizations will obey different equations but, at the 
same time, there will be differences in comparison with the vector case. The two polarizations 
of the relic gravitons when projected along the directions of the photon propagation will lead 
to a quadrupole term. The © and <g> polarization of the tensor mode are illustrated in Fig. 
E] which should be compared with the shaded area of Fig. |2J When the © polarization 
propagates orthogonally to the magnetic field direction its evolution equations are given by: 



d T Z @ + (ikfM + e')2© - ^d T h Q 



(C^ + A^S 



(t) 



2(l-/i 2 )cos2^ 

d T T(B + (ikfi + e')Tq 

e' 
2(1 -ii 2 ) cos 2<p 

d T S e + (ikfj, + e')«S ( 



AS(1 



I* 



C [Ai(l+f)-2fjAZQj?-l 



-^(C 2 A? 



A3)Eg 



AS(l+A* 2 )-C a (A?(l+/* a )-2^A a V-D 

3 



,(t) 

J © i 



,(t) 



16 



/eCAiA 2 



Vi 



fi' an <p ^ (t ) 



/! COS 2(fi 

tc/2. When the 



(6.17) 

(6.18) 
(6.19) 



where it has been assumed that a = (3 = —n/2. When the © polarization propagates 
orthogonally to the magnetic field direction its evolution equations are given by: 



d T Z® + (ik/j + d)Z® - ~d T h® = -e'CAiA 3 E® 
d r % + (ikfi + e')% = -e , CA 1 A 3 E§ ) , 
drS® + (ik/j, + e')«Sg) = 



e7oCAl A 3 £2^^lZZ E « 



(6.20) 
(6.21) 
(6.22) 



sin 2(p /j, 

Equations fl6U7l) . ( 16T81) and ( 16TT91) can then be compared to Eqs. f l6\20l . ( l?T2Ij) and ( 1fT22l) 
recalling that, now the magnetic field is directed along y. The two tensor polarizations read 
efj = (aiCij — bibj) and efj = (ctibj + djbi). But since it has been assumed that k = z we 
shall also have a = (1, 0, 0) = x and b = (0, 1, 0) = y. The addition of a magnetic field 
either along a or along b (i.e. orthogonally to k) is illustrated in Fig. [3] for the two tensor 
polarizations. The polarization © spans the shaded area bounded by the full lines. The 
polarization © spans the shaded area bounded by the dashed lines. The effect of having an 
extra source of circular dichroism along the a (or along the b axis) will be felt by both tensor 
polarizations as quantitatively established in Eqs. (16.170 — (16.19p and in Eqs. ( 16. 200 -1 )5.220 . 

The last step is to compute the evolution equations by averaging the source functions 
over the directions of the magnetic field. 



• (t) 



Cy(ji,tp,k,T) 

£g(/i,<M,r) 



Qzr' pi r27T rn r27C 

10Q o / dv dtp' sin(3d(3 dajf'(fj,, v, <p, <p', a,/3),(6.23) 

1287T^ J-l Jo Jo Jo 



rW/ 



12^ /-i dV C diP> So Sln m to da7 ^ (/ "' "' *' ^ "' ^ ' (6 ' 24) 
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Cy(fj.,(p,k,T) = / dz/ / tfy/ / sm(3d(3 da^'in, v, ip, <p',a,/3),(6.25) 

1287r^ 7-i 7o 7o 7o 

C\/{n,ip,k,T) = / di/ / <V / sm[3d[3 daJ^'fa, v, <p, <p', a, /3).(6.26) 

12o7r^ 7-i 7o 7o 7o 

The result for the evolution equations of the tensor polarizations with averaged sources is 
given by: 

d r Z + {ikii + e')Z - l -d T h = ^[C 2 (7A? - 5/ 2 A 2 ) - 6CA1A3 + 2A 2 ]£W (6.27) 



e' 



«9 T T + (ity + e')T = --[C 2 (7A? - 5/ 2 A 2 ) - 6(^3 + 2A 2 ]£« (6.28) 

<9 r S + {ikn + e')S = (6.29) 

These results extend and partially correct the results derived in [25]. The correction has to 
do with the source term of Eq. (I6.29P which vanishes exactly unlike stated in [25] because 
of an error in the azimuthal integrations. 

7 Concluding remarks 

An arbitrarily oriented magnetic field has been incorporated in the Stokes matrix of the 
last electron-photon scattering. The transport equations for the scalar, vector and tensor 
components of the brightness perturbations have been derived and studied in various physical 
situations. The obtained results pave the way for a consistent improvement of the available 
analytical and numerical tools used for the calculation of magnetized CMB anisotropies. The 
general treatment developed here is also expected to be relevant for the careful assessment 
of the level of circular polarization induced at last scattering by the presence of a small-scale 
component of the pre-decoupling magnetic field. 
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A Stokes and Mueller matrices 



The explicit form of the Stokes and Mueller matrices for arbitrary orientation of the magnetic 
field will now be reported. The four distinct entries of the Stokes matrix M(Jl, Q' , a, /3) 
appearing in Eq. (12.141) can be written, in explicit terms, as 



M 11 (Q,Q',a,P) 
, CAi + A3 



CAi-A 3 



1 — /i 2 Vl — v 2 + \w cos (if — a) cos (if' — a) 



2 
+ sin 2/3 



cos 2/3 



[W cos (f — a) cos (if' — a) — J 1 — /i 2 \/l — v 2 
li\J\ — v 2 cos (if — a) + uyl — fi 2 cos (if' — a) > 
+C A ly uz/ sin (if — a) sin (if' — a) + if e (A 2 \ sin /3 fiVl — v 2 sin (if — a 

+ fiv cos f3 sin (if — if') > , 



—v\j\ — fi 2 sin (ip' — a) 

M 12 (f2, fl', a, /3) = fi sin (if 1 — a) cos (if — a) 



(a.i; 



A 3 + AiC 



< n sin (ip' — a) cos (if — a) cos 2/3 + y 1 — fi 2 sin (if' — a) sin 2/3 > 



-(Aifi sin (ip — a) cos (<// — a) 



-i/eCA 2 

M 21 (fi,ft',a,/3) 



/i cos /3 cos (y/ — <£>) + y 1 — /i 2 sin /3 cos (<£>' — a) 
CAi + A 3 



(A.2) 



- a/1 — v 2 sin 2/3 sin (if — a) 



+I (CA 1 + A 3 ) 
-^(CA 1 + A 3 ) 



sin (if + if 1 — 2a) — sin (if' — ip) 
sin (y> + ip' — 2a) — sin (ip' — ip) 



cos 2/3 



ifMC 



' cos /3 cos (if' — ip) + Vl — v 2 sin /3 cos (if — 



oc, 



(A.3) 



M 22 (fi,ft',a,/3) 
4 



CAi + A 3 
4 



sin 2a(l — cos 2/3) sin (if' + ip) — A 3 cos (<£>' — 99) 



(1 + cos 2/3) [cos (<f' — ip) — cos 2a cos (if' + <£>)] 



CAi + A 3 



cos (ip' — if) + cos 2a cos (if' + if) 



-if e (A 2 cos /3 sin (if' — <p) . 



(A.4) 



Equations ( 1A.1J) . ( 1A.2J) . (1A.3|) and ( \AA\i have various interesting limits which depend upon 
the specific orientation of the magnetic field intensity. In the absence of magnetic field we 
have that 

A x -H, A 2 ->1, A 3 -+l, C^-l, /e-»0. (A.5) 
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In the limit denned by Eq. (|A.5|) . Eqs. f lA.ip - (IA.4p reduce to 



Mn(fi,(p,v : (p') = -yl - y?\J\ - v 2 - / uz/cos((/? / - (p), 

M 12 (/j,,(p,jy,ip') = fj, sin ((p' - if), M 21 (n,p,v,Lp') = -v sin {iff -tp), 

M 22 (/i, <p,v,(f/) = - cos (if' -ip), (A. 6) 

where, as already pointed out in section 121 /•* = cos?? and v = cos^'. The integrands 
appearing in the source terms of Eqs. (12.171) . (I2.18p . (12.191 and (I2.20p are given by 

7>(0, n', a, (3) = T n (n, Of, a, (3) + T IQ (n, M, a, P)Q(tf) 

+r w (n, n', a, p)u(n') + r IV (n, n', a, p)v{n'), (a.7) 

T Q (Q, Of, a, (3) = T QI (Q, Of, a, 0) I(fi') + T QQ (Q, Of, a, (3)Q(Q') 

+r QU (n, of, a, p)u(tf) + T QV (si, n', a , /3)v(n'), (a.s) 

JV(n, jy, «, /?) = r m (n, n', a, /3) i(fi') + r UQ (n, n', «, /3)Q(fi') 

+Tuu(n, Q', a, f3)U(Q') + T uv (^, fi', a, /3M1T), (A.9) 

jy(ft, jy, a, /?) = T V i(n, n', a, p)i{rt) + T VQ (n, n', «, /?)Q(o') 

+7W(^, n', a, (3)U(Q r ) + 7W(^, tt, a, f3)V(Q'), (A.10) 

where the matrix elements T%j are computed in terms of Eqs. f lA.ip . ( 1A.2I) . ( 1A.3[) and ( 1A.4J1 : 



Tn = 2[|M n | 2 + |M 12 | 2 + \M 21 \ 2 + \M 22 \% (A.ll) 

Tiq = 2[|M n | 2 - |M 12 | 2 + |M 21 | 2 - |M 22 | 2 ], (A.12) 

Tw = 2[M n M* 2 + AfiaKx + M 2l M* 22 + M 22 M* 2l ] , (A. 13) 

Tiv = 2[i(M 12 M* 1 - M* 2 M n + M 22 M* 2l - M 21 M* 22 ], (A.14) 

Tqi = 2[|M n | 2 + |M 12 | 2 - |M 21 | 2 - |M 22 | 2 ], (A.15) 

T QQ = 2[\M n \ 2 - \M 12 \ 2 - \M 21 \ 2 + \M 22 \\ (A.16) 

T QU = 2[M U M* 2 + M 12 M* ± - M 2l M* 22 - M 22 M* x l (A.17) 

T QV = 2[i(M 12 M* 1 - M* 2 M n - M 22 M* 21 + M 21 M* 22 ], (A.18) 

Tui = 2[M 11 M* 1 + M l2 M* 22 + M 21 M* X + M 22 M* 2 ], (A.19) 

T[/ Q = 2[M 11 M* 1 - M 12 M; 2 + M 2l M* xl - M 22 M{ 2 \ (A.20) 

Tuu = 2 [M n M; 2 + M l2 M* 21 + MaiM^ + M 22 M*,}, (A.21) 

7W = 2[z(M 12 M* 1 - M^jMii + Af^Mfi ~ ^1*2^21], (A.22) 

7^/ = 2[i{M 11 M* 1 + M l2 M; 2 - M 2l M* u - M 22 M* 2 ] , (A.23) 

T VQ = 2[i(M 11 M* 1 - M 12 M; 2 - M 2X M{ X + M 22 M* 2 ] , (A.24) 

T vu = 2[i(M u M; 2 + M 12 M* 2l - M 2l M* 2 - M 22 M* ± ], (A.25) 

Tw = 2[M* 2 M n - M l2 M* 2l + M 22 M{ X - M* 2 M 21 }. (A.26) 

In Eqs. flA.lip - flA.26P the explicit dependence upon the six angles has been suppressed only 
for sake of simplicity. The components of the incident electric fields in the local frame ei, e 2 
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and £3 can be related to the components off the electric field in the three Cartesian directions 

as 



£1 
E 2 



cos a cos (3E' X + sin a cos (3E' y — sin (3E' Z , 

— sin aE', + cos aE' 

cos a sin f3E' x + sin a sin (3E' y + cos /3£v 



(A.27) 



The incident electric fields E' x , E' and E' z can be related, in turn, to their polar components 

as: 



E' 



cos •&' cos (p'E'tf — sin (p'E' 



<p 
cos •&' sin ^'-E^ + cos (p'E' 

3/ 77i' 



¥" 



-sint?%, 



(A.28) 



where, as already spelled out in section[2]the direction of propagation of the incident radiation 
h' coincides with r' and (E'#, E' ) are the components of the incident electric field in the 
spherical basis. The relations between the outgoing and the ingoing electric fields is given 
by 



E4$,v,$',v',<y,P) 



r 



M 11 E#(#', ( p') + M 12 E 4 {#', ( p' 
M 21 E^' : (p') + M 22 E^' : (p' 



(A.29) 



where, M {j = A%(#, ip, ■&', if/, a, (3) are given by Eqs. flATTT) (TAT3]) . 



B Angular integrations 

The angular integrations of the collision term can always be performed either directly or 
with the help of the appropriate Rayleigh expansion. The results are however rather lengthy. 
Furthermore, in the case of the vector and of the tensor modes (i.e. sections Eland E]), the 
integrations over ip' change depending upon the the polarization of the vector (or of the 
tensor) wave. In the scalar case the results of the angular integration over ip' are given 
hereunder following the notation of Eq. (14. 6ft : 



T 
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2™ 2 (3/, 2 - 1) + fl 



(1 + /i 2 ) - 2(1 + /i 2 ) sin 2 (3 cos 2 (<p - a) 



+ 2vr 3-/i 2 + / ( 



;i+^ 2 



+ /lyl — /i 2 cos (ip — a) sin 2/3 

— 2 sin 2 /3(1 — /i 2 ) cos 2 (ip — a) + /iJl — fi 2 cos (<p — a) sin 2/3 



(B.l) 



r 



7T 



IQ 



[V 



1) 4(3/i 2 -l) + / ( 



4(1 - fx 2 ) sin 2 (ip - a) 



29 



+ (1 - /i 2 ) [cos 2(y? - a - /3) + cos 2(y? - a + /?)]+ 2 cos 2/3(1 + 3/i 2 






2/iy 1 — /t 2 ( sin (if — a — 2/3) — sin (ip — a — 2/3) 

0, 

87r/ e z4(l + L 1 " 2 ) cos/3 + /xy 1 — /t 2 cos (</? — a) sin/3], 



(B.2) 

(B.3) 
(B.4) 



27TZ/ 2 3(/i 2 - 1) + / c 2 



(/x 2 - 1) + 2 sin 2 /3(1 + (1 - fi 2 ) cos 2 (y? - a)) 



+ /iyl — /U 2 cos {ip — a) sin 2/3 



H-27rUl-yU 2 



+ /e 2 



(// 2 -l) + 2sin 2 /3 



sin 2 (</? — a) — /i 2 cos 2 (</? — a) 



+ /iyl — /t 2 cos (<f — a) sin 2/3 > 



(B.5) 



r, 



7T 



QQ 



1-z/ 2 ) 12(l-/i 2 ) + / ( 



2(/, 2 - 1) 



Tui 



Tuq 

Tuu 
Tuv 

Tvi 

Tvq 
Tvu 
Tvv 



- 2(1 + jj 2 ) cos 2(if- a) + (1 + /i 2 )(cos2(y?-a-/3) + cos2(y? - a + /3)) 

+ 6 cos 2/3(1 - /i 2 ) + 2fi\Jl - // 2 (sin (y? - a - 2/3) - sin (y? - a + 2/3)) j. (B.6) 

= 0, (B.7) 

= 8irf e is[(n 2 — 1) cos/3 + /tyl — /i 2 cos (y? — a) sin/3], (B.8) 

= — 27r/ 2 sin ((/? — a) < z/ 2 



+ 



4/i cos (ip — a) sin 2 /3 + \J 1 — /t 2 sin 2/3 
1 — fi 2 sin 2/3 — 4/x cos (<£> — a) sin 2 /3 
Anf 2 (i> 2 — 1) sin /3 sin (</? — a) 



1 — /i 2 cos /3 + 2/i cos (v? — a) sin /3 



0, 

—8insf c \Jl — /i 2 sin (</? — a) sin /3 

4tt/ c (z/ 2 



2/x cos /3 — a/1 — /i 2 cos (</? — a) sin /3 



+ 



2/i cos /3 + 3y 1 — /i 2 cos (y? — a) sin /3 

2tt/ c (z/ 2 - 1 
0, 



4/x cos /3 — 2\Jl — ji 2 cos (f — a) sin /3 



8"7TZ/< /i + cos /3/ 2 



/i cos /3 + y 1 — /i 2 cos (y? — a) sin /3 



(B.9) 

(B.10) 

(B.H) 
(B.12) 



(B.13) 

(B.14) 
(B.15) 
(B.16) 
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